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ABSTRACT 

The He I multiplet at 10830 A is formed in the incomplete Paschen-Back regime for typical 
conditions found in solar and stellar atmospheres. The positions and strengths of the various 
components that form the Zeeman structure of this multiplet in the Paschen-Back regime are 
approximated here by polynomials. The fitting errors are smaller than ~10~ 2 mA in the com- 
ponent positions and ^10 -3 in the relative strengths. The approximant polynomials allow for 
a very fast implementation of the incomplete Paschen-Back regime in numerical codes for the 
synthesis and inversion of polarization profiles in this important multiplet. 

Subject headings: line: profiles - Sun: atmosphere - Sun: magnetic fields - Sun: chromosphere 



1. Introduction 

One of the most useful near-infrared spectral 
regions for spectropolarimetric investigations of 
solar and stellar magnetic fields is that located 
around 10830 A, which contains both the photo- 
spheric line of Si I at 10827 A and the chromo- 
spheric lines that result from the He I 10830 A 
multiplet (Harvey & Hall 1971; Riiedi et al. 1995; 
Trujillo Bueno et al. 2002). Of particular interest 
is the fact that the 'blue' and 'red' components 
of this multiplet are sensitive to both the Hanlc 
and Zeeman effects, which offers a suitable diag- 
nostic window for investigating magnetic fields in 
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a variety of solar plasma structures, such as promi- 
nences (Trujillo Bueno et al. 2002), emerging flux 
regions (Solanki et al. 2003) and solar chromo- 
spheric spicules (Trujillo Bueno et al. 2005; Socas- 
Navarro & Elmore 2005). The He I 10830 A mul- 
tiplet originates between a lower term (2 3 Si) and 
an upper term (2 3 P2,i,o)- Therefore, it comprises 
three spectral lines (see, e.g., Radzig & Smirnov 
1985): a 'blue' component at 10829.09 A with 
Ji = 1 and J u = (hereafter Transition 1, or 
Trl for abbreviation), and two 'red' components 
at 10830.25 A with J u = 1 (hereafter, Tr2) and at 
10830.34 A with J u = 2 (hereafter, Tr3) which ap- 
pear blended at solar atmospheric temperatures. 

In a previous paper (Socas-Navarro et al. 2004, 
hereafter paper-I), we demonstrated that the de- 
termination of the magnetic field vector via the 
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analysis of the observed polarization in the He I 
10830 A multiplet must be carried out considering 
the wavelength positions and the strengths of the 
Zeeman components in the incomplete Paschen- 
Back effect regime. The theory of the Paschen- 
Back effect in this regime is detailed in Section 
3.4 of the book by Landi Degl'Innocenti & Lan- 
dolfi (2004) . The problem of finding the strengths 
and the splittings of the various magnetic compo- 
nents, arising in the transition from a lower and 
an upper term, both described in the L-S coupling 
scheme ( 3 S and 3 P, respectively, in the case of the 
A 10830 line), is reduced to the numerical diago- 
nalization of a number of matrices. The expres- 
sions of the matrix elements are given, in terms of 
Wigncr-symbols, in Eq. (3.57) of that book, or, in 
analytical form, in Eqs. (3.61a,b). As pointed out 
in paper-I, our results have been obtained through 
the adaptation to the fine-structure case of a com- 
puter program previously developed to handle the 
similar problem for hyperfme-structured multi- 
plets (Landi Degl'Innocenti 1978). 

Fig 1 shows the transition components with 
the linear Zeeman splitting (LZS) approximation 
and the rigorous calculation considering the in- 
complete Paschen-Back splitting effects (IPBS). 
The figure has been reproduced from paper-I, but 
corrects a plotting error in one of the Tr3 w- 
components. 

Significant differences between LZS and IPBS 
are to be expected for typical magnetic strengths 
in sunspots, simply because for fields stronger than 
about 400 gauss the Zeeman splittings of the J = 2 
and J = 1 levels of the upper term become com- 
parable to their energy separation and the pertur- 
bation theory of the familiar Zeeman effect is no 
longer valid. However, the discrepancies between 
the emergent Stokes profiles calculated assuming 
LZS or IPBS turn out to be significant even for 
field strengths substantially weaker than the level- 
crossing field of ~400 G (see paper-I). 

In this paper we provide polynomial approxi- 
mants that allow researchers to calculate Zeeman 
components and strengths of this interesting mul- 
tiplet in the IPBS regime. This will facilitate 
the generalization of numerical codes developed 
with the LZS approximation, such as the Milnc- 
Eddington code of Lagg et al. 2004, and to im- 
prove the efficiency of the IPBS code of Socas- 
Navarro et al (2004). Moreover, it may also be 



useful for solving the full non-LTE polarization 
transfer problem in the He I 10830 A multiplet 
within the framework of the incomplete Paschen- 
Back effect theory. 

Figure 2 demonstrates that our polynomial ap- 
proach to the IPBS regime in the He I 10830 A 
multiplet is sufficiently accurate for the calculation 
of the emergent Stokes profiles from a magnetized 
stellar atmosphere. 

2. Polynomial fits 

Let us characterize the deviation of the Zeeman 
pattern from the LZS as polynomials. For each 
Zeeman component, we have: 

AA/pb = AX L zs + V P {x) 

flPB = fLZS +Vf(x) , (1) 

where AA denotes the relative position of a given 
component (in mA), / represents the relative 
strengths which are dimcnsionless quantities (see, 
e.g., Eq [3.64] in Landi Degl'Innocenti & Landolfi 
2004), x is the magnetic field intensity expressed 
in kG, and y(x) are the approximant polynomi- 
als. Since the expressions for the splittings and 
the strengths in the IPBS regime have to con- 
verge to those of the LZS effect under the limit of 
weak fields, and since, under the same limit, the 
splittings are linear in the magnetic field and the 
strengths tend to a constant, it has to be expected 
that the differences y p and yf are of the form: 

y p = c 2 x 2 + C3X 3 + C42; 4 + c 5 x 5 
yf = d\x + d2X 2 + d%x 3 + d^x 4 , (2) 

with < x < 3. 

The following tables list the coefficients for the 
corrections y p (x) and yf(x) that need to be ap- 
plied to the LZS approximation to account for 
IPB effects. The last column of each table lists 
the maximum fitting error in the range considered 
(0 < x < 3). The components are sorted in the 
order of increasing wavelengths in the LZS regime. 

3. Conclusions 

The polynomial coefficients provided in this 
work can be used to generalize existing Zeeman 
transfer codes to treat the He I 10830 A multi- 
plet in the IPBS regime. As we pointed out in a 
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Fig. 1. — Positions and strenghts of the Zeeman components as a function of the magnetic field. In all 
panels the solid lines represent calculations considering IPBS while the dashed lines represent those using 
the LZS approximation. Note that, in the general IPBS case, the Zeeman components exhibit asymmetric 
displacements and strengths. Also notice that the Zeeman IPBS components have strengths that depend 
on the magnetic field. The central 7r-component of Tr2 has zero strength in LZS and is not plotted in the 
figure. The ordinate axes are given in Dopplcr units (corresponding to 105 mA) to make the figure directly 
comparable with that in paper-I. 



Table 1 
Component positions. TrI 



Transition C2 c 3 C4 C5 Max error 



1.685 3.072xl(T 2 -1.366xl(T 2 1.568xlCT 3 4.508xl(T 3 
1.685 3.073xl(T 2 -1.366xl(T 2 1.567xlCT 3 4.505xl(T 3 
1.685 3.073xl(T 2 -1.366xl(T 2 1.567xlCT 3 4.505xl(T 3 
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Fig. 2. — Stokes I and V profiles calculated according to a Milne-Eddington atmosphere for a 500 G field. 
Solid line: Strict IPBS calculation. Thin dashed line: LZS approximation. Thick dashed line: LZS corrected 
with the coefficients provided in this paper. The error induced by this approximation on the synthetic profiles 
is smaller than ^10~ 4 and is hardly noticeable in the figure. 



Table 2 
Component strengths. TrI 



Transition d\ g?2 c?3 ck Max error 

cr+ -1.040xl0" 2 2.326xl0- 4 1.577xl0" 5 -8.117xl0- 7 6.836xl0~ 7 

7T -1.280xl0- 7 -4.635xl0~ 4 2.127xl0- 8 1.606xl0- 6 1.034xl0- 7 

a- 1.040xl0- 2 2.331xl0- 4 -1.740xl0" 5 -4.781xlQ- 7 8.643xl0~ 7 
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Fig. 3. — Difference between IPBS and LZS component. Left panels: Position shifts (AXjpb — AXlzs)- 
Right panels: Relative strengths (fiPB — fhzs)- Top to bottom: Trl, Tr2 and Tr3. In all panels, the solid 
lines represent 7r-components and the dashed lines represent er-components (diamonds: cr + , squares: cr~). 
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Component positions. Tr2 
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Table 4 
Component strengths. Tr2 
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Table 5 
Component positions. Tr3 
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Table 6 
Component strengths. Tr3 
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Fig. 4. — Approximation errors. Left panels: Position shifts [(AXjpb — AXlzs) — Up]- Right panels: Relative 
strengths {{Jipb — Ilzs) — Vf]- Top to bottom: Trl, Tr2 and Tr3. In all panels, the solid lines represent 
7T-components and the dashed lines represent er-components. 
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previous work (Socas-Navarro et al 2004), neglect- 
ing IPBS results in significant errors in the cal- 
culation of its polarization profiles, even for field 
strengths much lower than the level-crossing field. 
Computer codes developed using these polynomi- 
als will be virtually as fast as those using the LZS 
approximation but the polarization profiles will be 
much more accurate. The actual improvement in 
CPU time obviously depends on the particular de- 
tails of the calculation. In our codes, the use of 
the polynomials is approximately a factor 20 faster 
than the full IPBS implementation. 

This research has been partly funded by 
the Ministerio de Education y Ciencia through 
project AYA2004-05792 and by the European 
Solar Magnetism Network (contract HPRN-CT- 
2002-00313). 
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